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Abstract 

We determine the widths of three-body decays of sleptons, l~ —>■ i^£^i^,£^iju,i~qq, in the 
presence of arbitrary slepton flavor violation and left-right mixing. These decays are important 
in scenarios in which the lightest supersymmetric particle is the gravitino, a generic possibility in 
models with gauge- and gravity-mediated supersymmetry breaking. Three-body decays have been 
discussed previously assuming flavor conservation and left-right mixing in only the stau sector. 
Flavor violation and general left-right mixing open up many new decay channels, which provide 
new avenues for precision mass measurements and may play an essential role in solving the standard 
model flavor problem. We present results for toy models with two-generation mixing, and discuss 
the implementation of these results in SPICE, a program that simplifies collider event simulations 
of flavor-violating supersymmetric models. 
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I. INTRODUCTION 



Fermion masses are one of the least understood parts of the standard model (SM). Even 
the charged fermion masses span over five orders of magnitude from the top quark to the 
electron, begging for a theoretical explanation. The majority of these masses and mixing an- 
gles are precisely constrained by experiment. Still, this wealth of data does not conclusively 
single out any theory of flavor. 

In the near future, this may change dramatically with the discovery of new particles at 
the weak scale. New particles may only deepen the mystery, as would be the case if a fourth 
generation were discovered. On the other hand, the masses and mixings of the new particles 
may be governed by the same principles that determine the SM fermion masses. In this case, 
rather than extending the fermion sector, the new particles will shed light on the existing 
fermion spectrum. 

Weak-scale supersymmetry provides examples of both possibilities. In pure gauge- 
mediated models, for example, squark and slepton masses are set by flavor-blind contribu- 
tions, with no connection to the SM fermion masses. However, in gravity-mediated models 
and hybrid models with both gauge- and gravity-mediated contributions, squark and slepton 
masses may receive contributions that are governed by flavor symmetries that also deter- 
mine the SM fermion masses [1-3]. The latter possibility leads to non-trivial flavor effects 
in high-energy experiments [4-6], opening the possibility for real progress on the SM flavor 
problem at the Large Hadron Collider, as has been emphasized recently by many authors 
(see, for example, Refs. [7-16]). 

Here we study the implications of flavor violation for the three-body decays of charged 
sleptons. Such processes are most relevant for colliders in models with a gravitino light- 
est supersymmctric particle (LSP) and a slepton next-to-lightest supersymmetric particle 
(NLSP), a generic possibility in models with both gauge- mediated supersymmetry break- 
ing [17, 18] and gravity-mediated supersymmetry breaking [19-23]. We consider a situation 
where the lightest particles are a gravitino LSP, several light sleptons, and the lightest neu- 
tralino, with the masses of each light slepton less than the mass of the lightest ncutralino. 
This results in the typical two-body decays of these light sleptons being highly suppressed 
or kinematically inaccessible, leaving the three-body decays as the dominant decay modes. 
In these cases, the three-body decays are also often the last visible step in cascade decays 
of squarks and gluinos, and so they impact nearly all supersymmetry searches and studies. 

Three-body slepton decays have been studied previously in an important, flavor- 
conserving case, where the authors considered cr erfi and jiR /^tti, with the fi 
a mixture of left- and right-handed staus [24] . These decays are characterised by two dis- 
tinct channels: a "charge-preserving" channel i~ —>■ with opposite-sign leptons and 
a "charge-flipping" channel i~ — > i'^i^i^ with same-sign leptons. Flavor-conserving three- 
body decays of squarks have also been considered [25], as have flavor- conserving three-body 
decays with sneutrinos as parent or daughter particles [26] . Our work generalizes the charged 
slepton analysis to the case of arbitrary lepton flavor violation (LFV) and arbitrary left-right 
mixing. In the presence of general LFV, any three-body decay of ii £j has up to 9 possible 
charge-preserving modes i~ — * ^J^k^m 6 possible charge-flipping modes , 
where ik,^m = r. In addition, LFV and left-right mixing bring additional compli- 
cations that are absent in the flavor-conserving case, including new processes mediated by 
Higgs and Z bosons, new final states with neutrinos and quarks, and new interference effects 
in charge-flipping processes. 
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These complications are well worth confronting, however, as there is a wealth of infor- 
mation in these branching ratios, which may shed light on the SM flavor problem. These 
decays, if present, are also relevant more broadly, for example, for supersymmetric searches 
and precision mass measurements. For example, in the flavor-conserving examples studied 
previously [24], the flnal states necessarily contain r leptons. Since these decay with miss- 
ing energy, they degrade searches based on energetic leptons, and they greatly reduce the 
prospects for precision mass measurements. With LFV, however, even if the lightest slepton 
is a stau, there may also be decay modes with two electrons, two muons, or an electron 
and a muon. Even if these branching ratios are suppressed, they may be the more obvious 
signals at colliders, and they may also provide better opportunities for precision mass mea- 
surements. It is therefore of interest to know the size of these branching ratios, and what 
determines them. 

In the following three sections, we begin with a general discussion of three-body decays 
and move gradually to more speciflc scenarios and concrete calculations. In Sec. II we 
present the new flnal states and new Feynman diagrams that are relevant to three-body 
decays once general flavor and left-right mixing are introduced. In Sec. Ill we discuss these 
results in more detail and explain the relative phenomenological importance of the various 
contributing diagrams in particular scenarios. In Sec. IV we then show concrete results 
in two toy models with two-generation slepton mixing to illustrate our results. Finally, in 
Sec. V, we present our conclusions and explain how our results have been incorporated into 
SPICE, a publicly available computation package for generating supersymmetric spectra and 
branching ratios in scenarios with arbitrary slepton mixing. 

We stress that, although we strive to give readers a intuitive feel for our results by 
considering concrete cases in the body of the paper, our analysis is valid for fully general 
LFV and left-right mixing. The complete, model-independent calculation is lengthy, but 
the full expressions for all three-body decay modes are given in a series of appendices. Our 
conventions and notations are deflned in Appendix A. These are consistent with those of 
Ref . [27] , where full details may be found. 

II. FLAVOR VIOLATION IN THREE-BODY SLEPTON DECAYS 

As mentioned in the Introduction, an important special case of three-body slepton decays 
has been discussed previously by Ambrosanio, Kribs, and Martin [24]. Motivated by pure 
gauge-mediated models, these authors considered the flavor-conserving decays cr — > erfi 
and jlR — >■ UTTi. These decays are mediated solely by neutralinos. The charge-preserving 
modes i~ — > iji'k'^rn shown in Fig. 1, and the charge- flipping modes —>■ t^^'^^^ are 
shown in Fig. 2. The charge-flipping mode is made possible by the Majorana nature of the 
neutralino, and detection of the charge-flipping mode would provide strong evidence that 
neutralinos are Majorana fermions. 

In the presence of more general slepton mixing, the gauge eigenstates e^, jiR, tr, e^, jli-, 
and tl mix to form six mass eigenstates, ii, i = 1, . . . , 6, with increasing mass, and the lepton- 
slepton-neutralino interactions are no longer flavor-diagonal. The neutralino diagrams of 
Figs. 1 and 2 are then modifled by the inclusion of 6 x 6 mixing matrix factors at the 
interaction vertices. 

In addition, new diagrams contribute. If the initial and flnal state charged sleptons 
contain left-handed components, there is the chargino-mediated decay to neutrinos shown 
in Fig. 1. This decay requires neither flavor violation nor left-right mixing. There are also 
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FIG. 1: Feynman diagrams for charge-preserving decays £^ ^j^k^m mediated by neutralinos 
and i~ — > iJi^kT^m mediated by charginos. 




FIG. 2: Feynman diagrams for the charge-flipping decays £^ ^j^k mediated by neutrahnos. 



charge-preserving decays mediated by neutral Higgs bosons and the Z boson, as shown in 
Fig. 3. The Higgs bosons mediate decays to same flavor and qq pairs, and the Z 

diagram mediates decays to same flavor uu, and qq pairs. 

Both the Ha and Z diagrams are present independent of LEV, but both require left-right 
mixing. In the case of the Higgs-rnediated diagrams, the Higgs couplings to sleptons have the 
form Hati^d-R, Hai}^iL, and Had-^^R- The first is generated by left-right mixing terms in the 
Lagrangian, and thus clearly requires left-right mixing to be non-zero. The other two come 
from D terms, and in the absence of left-right mixing, they are flavor-diagonal in the separate 
LL and RR sectors, even with the presence of LEV in one or both. Thus, in the absence 
of left-right mixing, the Higgs coupling to sleptons is flavor-diagonal in the 6-dimensional 
space of left and right sleptons and does not contribute to three-body decays. Similarly, the 




4-9 (■k^Vk,q 



FIG. 3: Feynman diagrams for charge-preserving decays mediated by Higgs and Z bosons. The 
Higgs scalar Ha may be any of the neutral Higgs bosons: , H^, A^, or (in Feynman gauge) G^. 
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Z couples to the LL and RR sleptons separately, so in the absence of left-right mixing, its 
couplings are also flavor-diagonal. Note, however, there are no photon-mediated diagrams at 
tree level, as the photon couples with equal strength to the left- and right-handed sleptons, 
and so its couplings are always flavor-diagonal. Also note that in the case of non-zero left- 
right mixing but no flavor mixing, the and Z diagrams exist, but their contributions 
are only appreciable for certain mass orderings. Namely, only when neither of the two mass 
eigenstates produced by left-right mixing posses kinematically allowed two-body decays will 
the three-body Ha and Z contributions be relevant. 

The charge-flipping decay is unaffected by the Higgs and Z boson diagrams, but in the 
presence of flavor violation, the fermions in the final state may be interchanged. This implies 
a new interference effect in the charge-flipping mode which is absent in the flavor- conserving 
case. In principle, there should also be a similar interference term in the neutrino decay 
modes if neutrinos possess Majorana masses, but such a term is negligible in the limit of 
zero neutrino mass. 

In summary, for charge-preserving decays to charged leptons of different generations, the 
analysis of Ref. [24] may be straightforwardly modified to the general case. The only modifi- 
cation needed is the insertion of rotation matrices from flavor eigenstates to gauge interaction 
eigenstates. For charge-preserving decays to charged leptons of the same generation, how- 
ever, the Higgs and Z diagrams introduce new contributions to the decay width. There are 
also new charge-preserving decays to neutrinos, mediated by charginos and Z bosons, and to 
quarks, mediated by Higgs and Z bosons. Finally, in the case of the charge-flipping decays 
to like-sign leptons, the daughter leptons may be interchanged and thus the decay width 
must include an interference term between these two diagrams; however, there are no Higgs 
and Z contributions to this mode and so there are no further interference terms. 



III. ANALYTIC RESULTS AND QUALITATIVE OBSERVATIONS 

Given an understanding of the qualitatively new features introduced by flavor and left- 
right mixing described in Sec. II, we can now calculate the resulting decay widths. The 
full results arc lengthy and arc given in the appendices. In this section, we summarize the 
results of our calculations and provide qualitative insight into the relative importance of 
each contributing mode. 

Our slepton interaction Lagrangian is 



/'int 



+aiftj,h'' + affi*i,H' + taffiU.A' + ta^i*e,G' , (1) 

where the coefficients j3, 7, (, and a contain gauge and Yukawa couplings along with the 
mixing matrix elements for sleptons, Higgs bosons, neutralinos, and charginos. These coef- 
ficients are defined in Appendix A. Here, the indices are summed over the six slepton, three 
lepton and neutrino, four neutralino, and two chargino mass eigenstates. We have kept only 
the terms relevant to the three-body decays in question. 
Our lepton interaction Lagrangian is 

(^) 

= ^4 (/iO sin Oh - cos Oh + i^°7^ sin (5 - iG^j^ cos (3) Ik 

+—^Z, (1 - 2 sin^ 9w) 47^4 ■ (2) 
2 cos 9w ^ ' 
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The amplitudes for charge-preserving decays to charged leptons 1^ — > £^ ^+ are 



M^o = -iu (p,) {PI*,Pr + PI*,Pl) ' [Pl^mPL + P%^Pr) V (Pm) (3) 

[Pk ~ Pi) ~ 

Mz = TT^l ^ (^^^) (^i + Pi) (2 sin' - Pl) V (pm) 6km (4) 

2 cos t^iy (pi — pj) — rriz v » ' v / 

- (2) {£) ■ n 1 

MhO = -7= 72 [pk) V {pm) hm (5) 

V ^ [Pi — Pj) — TTlfj^o 

■ (3) {() a 1 

Al/fO = 72 [Pk) V [Pm) dkm (6) 

V ^ \Pi ~ Pj) ~ n^ffo 

■ (4) {£) ■ o 1 

AIao = 72 (Pfe) 7 (Pm) dfem (7) 

■ (5) W Q -, 

^(T}/|/;^COS/? 1 _ . 5 / w /ON 

A^GO = 7^ 7 ^2 [Pk] 7 V [Pm) hm , (8) 

V2 (Pi - Pj) - m| 

where the indices i, j, k, and m correspond to the subscripts specifying the mass eigenstates 
of initial and final state particles in the decay mode i~ — >• ijij^i^, and 6km is the Kronecker 
delta function. Here and throughout the rest of this section we suppress these indices on 
matrix elements and decay widths. 
The overall decay width is given by 



+2 Re 



.6=1 a b=l b 



(10) 



= r^o^o + 6km i^ZZ + ^HH + '2^Zx° +'^^Hx° +'^^Zh) , (H) 

where /pg^ is the integral over three-body phase space discussed in Appendix B which in- 
cludes a sum over fermion spins, and we use an intuitively obvious notation for widths, so 
that, for example, ^x°^o is the partial width from J2ab-^x2.'^l°- 

To develop a qualitative understanding of which matrix elements are typically dominant 
and sub-dominant, we first note that there is a suppression to the Higgs and Z modes from 
left-right slepton mixing. In fact, the explicit Higgs and Z couplings to sleptons depend on 
the left-right slepton mixing squared. 



(a) A2) _ ^ . 2) 



In particular, for /i not significantly larger than the slepton mass, this is roughly the same 
order as Yukawa suppression. We thus expect that, at the matrix element level, the Higgs 
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or Z radiated modes are generically suppressed by two powers of left-right mixing relative 
to the neutralino and chargino modes. 

Second, to simplify phase space factors, we assume a typical amount of mass squared 
splitting l\w?^lrv?^ < 5%; we also neglect lepton masses when they are sub-dominant. Under 
these assumptions, the difference in propagator structure between the neutralino mode and 
that of the Higgs and Z modes becomes apparent. The neutralino propagator is inversely 
proportional to [pk ^ Vif ~ ""^lo ~ "n^f — ""^lo, while the Higgs (Z\ propagator is inversely 
proportional to [pi — pj)'^—m'jj^ ^ ^ m| —m'j —Tn?Ha zi which simply reduces to —fn'jj^ ^ if the 
slepton masses are not too far above rriz- Thus, in models where the lightest neutralino and 
light slepton masses are close, the lightest neutralino pole contribution will be enhanced over 
the Higgs and Z mass-suppressed contributions. Conversely, as the slepton and neutralino 
mass scale grows, the Higgs and Z contributions will drop off more slowly than the neutralino 
contribution. 

The three-body decays with two neutrinos or two quarks have much the same form 
as above with minor changes (refer to the appendices for details). These changes are as 
follows: For the neutrino mode, Pl^j^ —>■ jaik and P^^f, — > 0, and the lepton masses are 
set to zero, hence simplifying the phase space calculation and removing the off-shell Higgs 
contribution. For the quark modes, there is no off-shell fermion intermediary, and the Higgs 
and Z modes are only modified with adjusted couplings and quark masses. The quark 
modes are, however, enhanced by the color factor and a sum over light flavors. Because of 
these changes, however, all of these modes are suppressed relative to the di-lepton mode by 
left-right mixing. The quark modes only have contributions from the Higgs and Z diagrams, 
and thus are suppressed by left-right mixing as noted above. In the majority of models with 
a hght slepton NLSP, the in gauge eigenstates are generically hghter than the II gauge 
eigenstates;^ for such a model the neutrino modes arc suppressed since the neutrino only 
couples to El gauge eigenstates while the light sleptons are primarily iji states. Indeed, in 
such a model 



so the neutrino and quark decay widths are all suppressed by (rnj^j^/m'jj relative to the 
dilepton modes. 

Finally, for the charge-flipping di-lepton decay, the matrix element is given by 



M§ = -in ip,) {f3l*^Pn + PI^Pl) % _ ^° {PI%Pl + (3]:,,Pr) v (pj (15) 




(13) 






(16) 



where the negative sign comes from Fermi statistics. The decay width is then 




(17) 



with 




(18) 
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and Ckni is a phase space factor which has a value of 1 for two outgoing leptons of different 
flavor and 1/2 for two outgoing leptons of the same flavor due to indistinguishable particle 
statistics. Fn and T22 both have the same basic form as the result from [24], where the only 
change is the insertion of flavor- mixing coefficients and, for r22, interchange of £k and 
between the two terms; the new is presented in Appendix D 4. We find that this charge- 
flipping decay width is of the same order of magnitude as the charge-preserving di-lepton 
mode, though the flavor structure is markedly different: this is discussed in the next section. 



IV. TWO-SLEPTON MIXING: ILLUSTRATIVE EXAMPLES 

To validate our results and investigate their phenomenological implications, we examine 
some simple cases of two-slepton mixing. In these examples we consider spectra with fairly 
degenerate sleptons and lightest neutralino. The slepton-neutralino degeneracy is motivated 
by simple gauge-mediated supersymmetry breaking scenarios with not too many messengers. 
Larger splittings are, of course, possible if there are many messengers or in other frameworks, 
such as minimal supergravity. The sleptons are taken fairly degenerate so that large mixing 
angles are consistent with low energy constraints and so may be considered. We note, 
however, that these examples are merely illustrative, and our results are valid in any chosen 
framework with arbitrary mass splittings. 

First, we consider cr — jlu mixing parameterized by 

/ cos ^12 sin ^12 
\ —sin ^12 cos ^12 





(19) 



Fig. 4 shows the flavor- violating decay widths (.2 if£~£'^ with electron and/or muon 
leptons as a function of the mixing angle 9i2- Here both leptons are taken to be explicitly 
massless and left-right slepton mixing is set to zero, thereby cutting off the Higgs and Z 
modes, ii and £2 are assigned masses of 100 GeV and 105 GeV, respectively, and the lightest 
neutralino is given a mass of 110 GeV (heavier neutralino contributions are small in this 
case) . 

The left plot in Fig. 4, which shows the decay widths for the charge-preserving channel, 
demonstrates the typical structure of two-slepton mixing: all decay widths are at most tt 
periodic, the sleptons interchange roles (/x~ — > e~ and e"*" — >• /x+) at 6*12 = tt/2, the e~e+ and 
fj,~fi~^ modes are equal at all mixing angles, and at 9i2 — tt/A all decay widths are equal. 
In contrast, the right plot in Fig. 4, which shows the charge-flipping channel, demonstrates 
a different flavor structure. Since the e~^~ and /x~e~ modes in general contribute to the 
same decay width, all widths are 7r/2 periodic. Also, the e~e~ and modes are again 

equal at all mixing angles. In addition, note that the ii~e~ mode drops to zero at 9^2 = 7i'/4, 
because the charge- flipping mode has two diagrams at tree level which cancel at 9i2 — tt/A 
for the fjre~ mode. Then, as expected, the e~e~ and iJi,~iJr modes have decay widths that 
are half the total width at 612 = 7r/4. 

Next we consider cr — tr mixing, with massive leptons and a mixing angle ^13 given by 




cos 6'i3 sin ^^13 
— sin ^13 cos ^^13 




(20) 



Figure 5 shows the widths for £3 ifi where the leptons are electrons and/or taus. 
ii and £3 have masses of 100 GeV and 105 GeV, and lightest neutralino mass is again 110 
GeV. Left-right mixing is again neglected. 
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FIG. 4: Di-lepton decay widths as a function of mixing angle between smuon and selection with 
left-right mixing neglected. The left plot shows the charge-preserving channel, the right plot shows 
the charge-flipping channel. In both plots the leptons are taken to be massless, nip = 100 GeV, 



m.- = 105 GeV, and m^o = 110 GeV. 
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FIG. 5: Di-lepton decay widths as a function of mixing angle between stau and selectron with left- 
right mixing neglected.. The left plot shows the charge-preserving channel, the right plot shows 
the charge-flipping channel. Here the leptons are both massive, m^^ = 100 GeV, m^^ = 105 GeV, 
and m^o = 110 GeV. 



The first notable feature in Figure 5 is the separation between the e~e^ and r~r^ modes: 
the e~e^ mode is almost unchanged numerically from the previous case, while the r~r^ mode 
is suppressed by phase space constriction, as expected. Likewise, the r~e^ and e~r^ modes 
are also suppressed relative to the massless case, but certainly less suppressed than the r~r^ 
mode. Note that the total decay width is no longer constant because of the nonzero tau 
mass. 
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V. CONCLUSIONS AND IMPLEMENTATION IN SPICE 



We have determined the decay widths of three-body slepton decays, including both 
"charge-preserving" processes l~ — > i''i~£~^ and "charge- flipping" processes £~ t^l'l" , 
in the presence of arbitrary lepton flavor violation (LFV) and arbitrary left-right mixing. 
Such processes are particularly relevant in scenarios with a gravitino LSP and a slepton 
NLSP, where they are typically the dominant decay of some of the non-NLSP sleptons and 
are present as the last step in many SUSY cascade decays. 

Arbitrary 6x6 slepton mixing leads to many new diagrams and new decade modes, as well 
as new interference effects. Our results are fully general, but we have illustrated them for two 
simple cases with 2-generation LFV. The branching ratios to new modes may be significant, 
even accounting for low-energy constraints on LFV. The new modes may substitute final 
state electrons and muons for taus, potentially providing more obvious signals and better 
prospects for precision mass measurements. In addition, the measurement of these branching 
ratios may help pin down the supersymmetric fiavor parameters and shed light on the SM 
fiavor problem. We note that our illustrative examples have assumed fairly small slepton 
mass splittings so that large mixing angles are consistent with low energy constraints. It is 
also possible, of course, that the splittings are larger and the mixing angles are smaller. The 
calculations presented here are valid also in those cases, and it would also be very interesting 
to determine whether such small mixings may be established through collider studies. 

The full calculation presented here is lengthy, but has been incorporated into SPICE: Sim- 
ulation Package for Including Flavor in Collider Events, a publicly available code. SPICE 
takes as input a fiavor-conserving model, such as minimal gauge-mediated supersymmetry 
breaking, and arbitrary lepton fiavor- violating parameters, and generates the full supersym- 
metric spectrum, complete with fiavor-violating branching ratios. The three-body decays 
discussed here are included in SPICE. Details on obtaining and running SPICE are described 
in the SPICE user manual [27]. 
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APPENDIX A: LAGRANGIAN 



As given in Eq. (1), the relevant flavor- general interactions for this study are given by 
the Lagrangian 



pint 



(2) 



la. 



ij ^i^3' 



'ij ^1^3 



(2) 



id t 



(Al) 
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The coefficients are 



^ {goi. + 9' 01,:) uit - yPoiXll 



3^^^ = 

_ \-g cos Ol uS* + sin Ol ufll, 
[ g sm Ol UH + yl cos Ol UH^ ^ 



a = 1 
a = 2 



(^^(l-tan^MsK^H + Z?) 



(A2) 
(A3) 
(A4) 

(A5) 



ffm^^singjf \ ^(e)*^{e) 



\ 2 

+ (gm^ tan^ sin(^, + /5) - t/ig^t/S. 
y mw'Cosp J ' '■' 



+ g'mvK tan cos{Bh + P) ^ 

\ mw cos p 



^(,+4t^ff)(uS'u%-ui2ifjS 



2mw 



(A7) 
(A8) 
(A9) 



^i.c"* tfie slepton and sneutrino mixing matrices, respectively. Tfiese relate the 
slepton and sneutrino gauge eigenstates, 



(AlO) 
(All) 

to the mass eigenstates z/j through the relations ii = U-^^ia and z/j = Ui^^^Ua. The 
mass eigenstates are defined in order of increasing mass. The neutralino gauge eigenstates, 

-0° — (^—iB, —iW, tp% tp^^ , are related to the mass eigenstates Xi by Xi = Ojai'a- Similarly, 



the mixing of the negative charginos 



Xi \ / cos6'i sin6'L 
X2 J V - sin Ol cos Ol 



-iw- 



Finally, the neutral Higgs doublet is given by 



(A12) 



(A13) 
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where 

„ / cos Oh sin 6'ii- \ ,7-, / sin B cos B\ / a -, . \ 

ReH^\ -a -3 and = ■ n \ ■ (A14) 

^ I — sm 9h cos 9h I \ — cos p sm fj j ^ ' 

These coefficients are consistent with those of Ref . [27] , where a more complete discussion 
may be found. 



APPENDIX B: PHASE SPACE INTEGRATION 



In the three- body decays — > ijfkfm, we label the initial and final state particle masses 
rrii, rrij, m^, and rrim- We will maintain these indices throughout the appendix, with i 
corresponding to the parent slepton, j to the daughter slepton, k to the daughter fermion 
with the same sign as the parent slepton, and m to the other daughter fermion. The decay 
widths are the squared matrix elements integrated over 3-body phase space: 



\M\' ^ 



1 



dEj 

Jet Je, 



dEm\M\ 



rrii 



2567r3 



ry+ 

/ dx / dy\M\' 

Jx- Jv- 



where the dimensionless quantities 



2E, 



X — 



rrii 



and y 



2E„ 



rrii 



have limits 
x^ - 

y± - 



1 + ^fc - i^j + ^m) 

1 



2{l-x + rl) 
and we have defined 



(2 - x) (r^ -x)± Y^a;2 -4r^A5 (l - x + r^, r 



r 



rrii 



rrii 



rrik 

rrii 



rrir. 



rrii 



X{x,y,z) 



1 2 I 2 I 2 

2 I 2 I 2 
X +y + z 



2xy — 2xz — 2yz . 



(Bl) 



(B2) 



(B3) 
(B4) 

(B5) 



(B6) 

(B7) 
(B8) 



Throughout this work, we implicitly assume that the matrix element squared has been 
averaged and summed over initial and final state spins, respectively. 
The squared matrix element has the form 



(B9) 



a,b 



where Aia and Aib are matrix elements corresponding to Feynman diagrams with inter- 
mediate particle a and b. To evaluate the decay widths, we choose to integrate over y 
analytically. The decay widths are then written in terms of dimensionless integrals over x, 
which are performed numerically. 
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In integrating over y. it is convenient to note that all terms Aia-M-l have a numera- 
tor that is either constant or linear in y and a denominator that is proportional to the 

product of two propagators, each of which is either of the form (l — a; + — or 

(^x + y — 1 + Tj — r^^j , where ra,& = mafi/rrii. We therefore define the following dimen- 
sionless functions of x: 

r 

f{x,ra,rb) = / 

Jv 



y+ ydy 

y_ {l-x + rl){l-x + rl- rl) {1-x + rl- rl) 



' {x,ra,rb) 



/i(^) (x, r„, Tfe) 



k {x, r-b) 



Jv 



{l-x + rl-rD {l-x + rl) 
y+ dy 



y- {i-y + r^-rl){l-x + rl-rl) 
l{x,rl-l- rl) 



(l-x + rl-rl) 

y+ ydy 



{l-x + rl- ^b) 

y+ ydy 



ry+ dy 
V- {x + y — 1 -\- r'j — r'l^ 

I {x^x — 1 + r'j — j — I {x^x — 1 + r'j — rfj 



{rl - rl) 



y+ ydy 



y- (x + y — 1 + rj — r^j 

y+-y- 



{y+ + X - 1 + r| - j {y^ + x - 1 + r| - j 

Jv 



(BIO) 



(Bll) 



y- {^~y + rl,-rl){l-x + rl-rl) 
= (l + - rt) f''> {x, r„, r^) - k {x, r^) (B12) 

9^^Hx,ra,n) = ^U- - 

Jy- [x + y-l + r^ -rl){l-x + ri- r^) 

l(x,x-l + r] -rfj 



(B13) 



'y- (x + y - 1 + r'^ - rfj {1 - X + rl - r^) 
- k {x, rb) + (l-x- r| + r^) g^^^ {x, r^, r^) (B14) 

/iW (2;, r„, rb) = "Y 

Jy- ( a; + 1/ — 1 + r| — 1 ( a; + y — 1 + r| — 1 



(B15) 



'y- (^x + y — 1 + r| — r^^ (^x + y — 1 + r| — r^^ 
= / (x, X - 1 + r| - r^^) + (1 - X - r| + r^) h^^^ {x, ra, r^) (B16) 

.(1). V fy+ dy 
> {x, ra) = / 
Jy- 



(B17) 



y- (x + y-l + r]-rl) 
l{x,x-l + r]-rl) + {l-x-r] + rl)i^^^{x,ra) , (B18) 
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where 



k {x, Th) 

I {x, z) 



y+-y- 



l-x + rl-rl 



In 



y+ + z 



y- + z 



(B19) 
(B20) 



These functions will appear frequently in the decay widths to be discussed below. The 
first function / contains an extra factor in the denominator to conform to the notation used 
in Ref. [24]. The other functions are simply integrals over y with all relevant combinations 
of propagators in the denominator and a numerator either constant or linear in y. 

APPENDIX C: T (iT ^ £7^-^+) 
1. Total Width 

For the charge-preserving case, the matrix elements for all contributing modes are pre- 
sented in Eqs. (3)- (8). The total decay width is 



j k m 



1 



1 

~ 64^ 
+2 Re 



(CI) 



b=l a b=l b 



(C2) 
(C3) 



where Skm is the Kronecker delta function. For F^o^o the general case of independent k and 
m is taken; for the remaining widths, A; = m is assumed to simplify the expressions. 



2- XX Contribution 

The neutralino width corresponds to the width given in Ref. [24] with the addition of 
flavor- violating vertices. The decay width is given by 

r«-^i:Ec^|;«4?,.. (C4) 

t=l a,b=l 

where t labels the coefficients and integrals, and a and b label the neutralinos. The dimen- 
sionless integrals are defined by 

4gxg = £^ - 2r2) (l - X + rl) (r^ - x) f (x, r„, n) (C5) 

4?x° = ^"-^b j dxix- 2rQ {r^ -x) f (x, r„, n) (C6) 
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2r^r6 J ^ dx(x- 2rfj (l-x + rfj f (x, Ta, n) 
4?xg ^ ^^^'^^ / dx(l-x + r^) (i?^ - 2;) / (a;, r„, rt) 
'^rkTmran j dx(l-x + rlj f {x, Va, n) 
2rfcrm y (1 - X + ry / (x, ra, n) , 



r(5) 



7-(6) 



with coefficients 



Aa 



r>U n2 nl n2* , ^2* ^1 n2 nl* 
r'aikl-^ ajml-^bikl'^bjm ' 1-^ aikl-^ ajml-^hikrbj 

2 Re PaikPajmPbikPbjm ~^ Paikf^ajmPbikPbj 

— 2 Re PaikPajmPbikPbjm PaikPajmPbikPbj 

B^* B"^ B^* 

r'aikr^ajmr^bikr'bjm 

B^* B^ B"^ B"^* 
aikr^ aimr^bikr^bjm 



^(3) 

Aa Af) 



r<y^> — 

^<^o^o — 

Aa A5 



C 



2* 



(5) 
y" y" 

Aa A/) 



^(6) 
•-^vOyO 

Aa A5 



-4 Re 
-4 Re 



(C7) 
(C8) 
(C9) 
(CIO) 

(Cll) 

(C12) 
(C13) 
(C14) 
(C15) 
(C16) 



ZZ Contribution 



The Z boson contribution to the decay width is 



^7,7, — 



ITii 



5127r3 

where the dimensionless integrals are 



cos 6' 



w 



S^Mt) At) 

t=i 



^ZZ 



r(2) 
^ZZ 



+ 



^2 
^k 



dx (2a; - 2 - {x,rz) 
'.-{2 + rl)x-2r^^+3rl + rlri 



,•(1) 



{x, rz) 



dxi^^^ 



x,rz_ 



+ (x-3-r|) 



i^^) {x,rz) , 



and the coefficients are 



C^zz = 8 sin^ Ow-^ sin^ Ow + 1 
Cgl = 8 sin^ Ow-^ sin^ . 



(C17) 



(CIS) 

(C19) 



(C20) 
(C21) 



4. HH Contribution 



The purely Higgs-mediated contribution actually consists of several pieces with similar 
phase space structure. Furthermore, it simplifies into the sum of contributions from the 
real and pseudoscalar Higgs bosons, since the interference term between a real scalar and a 
pseudoscalar vanishes. 
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The width is given by 



1 



^HH — o {ChPh°^hOhO + C HO HO I HO + 2C/i0^^o/^ojj-o 

zobTT-^rrii 

+ Caoao^aoao ^G°G°Ig°G° '^^A°G°Iaog°) ' 

where the dimensionless integrals are 

Ikn. = £^ dx i^'^ (x, rnj + {x - - 2rl) i^'^ (x, r^J 

IkH, - dx h^'^ {x, rH„, rnj + {x - - 2rl) h^'^ {x, th^.th,] 

^kH, ^ / ^ ^^^^ ^Ha,rH,) + (x-l + r|) /i^^) (x, rH,,rH,) , 
where a ^ b, and the coefficients are 

dri ft 



ChOh° — 

Choho = 

ChOHO — —Re 

Caoao — 



(3) W a 
a)i'yl' cos Oh 



(4) (^) • Q 



(2) (3)* (<?)2 . ^ ^ 

c^jVc^jV Uk sm 9 H cos Oh 

2 



C^OGO = a\fyi^^ COS f3 



C 



AO GO 



-Re 



(4) (5)* (<?)2 . ^ ^ 

(^ji (^ji Vk sin/3cos/3 



(C22) 

(C23) 
(C24) 
(C25) 
(C26) 



(C27) 

(C28) 
(C29) 
(C30) 
(C31) 
(C32) 



5. Z-x^ Contribution 

The width from the interference term is 



a(2) 4 2 

rrij g(,ji ^ ^ ft) 
'BUn^cosOwtih ^^'^ 



where the integrals are 



^zlo = £^ dx[2(l-x + rl) ^(2) {x, rz, n) 

+ {2x - 2R^ + rl- r^rl^ g^^^ {x, rz, r^) 



rkn I dx g^'^^ {x, rz, rb) - {l + x - rfj g^^^ {x, rz, n) 

J X — 

g'-'^^ {x, rz, rb) + (l-x- r|) ^f^^^ {x, rz, n) 



7-(2) 



(3) /'^+ 

^zxl = ^k^^ I dx 



= ^^r^^^ [{3-2x + r])g^'^{x,rz,n) 



(C33) 



(C34) 
(C35) 

(C36) 

(C37) 
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and the coefficients are 



;(2 sin^ dw - l) PlikPlU + 2 sin^ Ow^LPI-^ 



ego - Re 



ego = Re 
ego = Re 



sm dw 

2 sin^ dwPLPl*m + (2 sin^ dw - l) /^L/S^; 
' " '/3L/36;;n + (2 sin^ Ow - 1) /^L/^S 



im 



2 sin^ ^w^/ 



6. Hx^ Contribution 

The width from the i?x° interference is 



where the dimensionless integrals are 



tH{1) 




1 dx 

'x- 




{x,rH^,n) + (x-R^-2rl)g^'^ 


{x,rH,,rb) 


jH{2) 




/ dx 

Jx- 


(2x 


-i?2-2r2)^« {x,rH^,n) 






^HaXl 


= n 


1 dx 

Ix- 




(x, TH^.n) + (x - 1 + rfj g^^^ (x, 






^HaXl 


= r-k 


I dx 

Jx- 




rfj g^^'> {x,rH,,rb) , 







and the coefficients are 



e 



(1) 

h<'x°t 
^(2) 



C 



(1) 



^(2) 



Re 
Re 
-Re 
-Re 
Re 
Re 
-Re 
-Re 



(2) (i) ■ n ( ol o2* I n2 nl* 
^ji Uk Sl'^ yPbikPbjm ' PbikPbjm 

(^fi'Vk^ sin Oh (Plikf^ljm + PlikPbjm 

cos Oh {filk^l^m + ^L^l^) 

'^ji^Vk ^ (^(^bikPbjm + Pbikf^bjm) 

'^ji^Vk ^ (^PbikPbjm ~ l^bikl^bjrn) 

^ji^yi "* sin P (PbikPbjm ~ f^bikf^bjm) 
^^iVk ^ COS/3 ij^lj^kPlj^ — Pbikf^bjm 
^ji^Vk ^ cos/5 [^f^likf^ljm ~ PbikPbjm 



7. ZH Contribution 

The ZH interference decay width is 



^ZH — 



cos 9w 



{Czh°Izh° + CzHoIzH° + CzAoIzAO + CzGoIzG°) ' 
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where 



and 



'ZHa 



^ZHa 



Tfc / dx h^'^^ {x, rz, thJ - xh^'^'^ {x, rz, r^J 
Tk f ^ dx (l- r|) h^^^ {x, rz, thJ , 

J X — 



CzhO 
CzHO 
CzAO 
CzGO 



-Re 
Re 
-Re 
Re 



afyi'^ sin Oh {ism^Ow-l) 

' (4) {(■) ■ Q 

(5) (i) n 



(C56) 
(C57) 

(CSS) 
(C59) 
(C60) 

(C61) 



APPENDIX D: T [ir ^ lfe^£-) 
1. Total Width 

For the charge-flipping decay, the calculation is complicated by the fact that the same-sign 
daughter leptons create an interference term. It is convenient to break the matrix element 
into two parts, as shown in Eqs. (14)-(16). It is correspondingly convenient to separate the 
decay width into three terms 



k ^m) — Ckm (ril + T22 - 2121) , 



where 



r = 



1 



j2 RelM'i^Mi'^' 



(Dl) 



(D2) 



and Ckm is 1 when the two outgoing leptons are of different generations and 1/2 when they 
are of the same generation. 



2. Til Width 

From inspection of the matrix element, Ai^a^ is identical to At^o except for the coefficients 
in the ajm vertex. Then Fn is identical to F^o^o in the charge-preserving case with the 
substitutions 



Pi 



xjm 



and with the identical substitutions for the complex conjugates, where x — a,b. 



3. T22 Width 

From further inspection of the matrix elements, M.^a^ and M.^^^ differ only in the inter- 
change of the two outgoing leptons. Thus F22 may be obtained from Fn with the interchange 
k m. 
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4. Tai Width 



The width from the interference term is given by 

1 



r2i = 



/ y Re 



2567r3 



EE Re 

t=l a,b=l 



where the integrals are 



t(1) 



j(2) 
XaXft 

j(3)_^ 

XaX^ 

t(4) 

XaX^ 
XaXfo 

7(7) 
t(8) 

AaA^ 



£^ { [2rlrl - X (l + r^) + i?^] /(^) (x, r„, r,,)} 

- (l - X + rQ (x,r«,rb) 

rkVm I dx {l + r| - - r^) /^^^ (x, r^, r^) 

J X — 

rmTb / rfa; \R^ ( r-fe) - /^^^ (a;,ra,rfo) 

Tkn j dx [/(^^ (X, Ta, Tft) - 2r^/(^^ (x, Ta, Tfe) 

?^m?^a y dx (x- 2r^) /(^) (x, r„, n) 
TkTa j dx (r^ - /^^^ (x, r„, n) 

rb) , 



and the coefficients are 
D 



(1) 

■0^0 



Aa Af, 

r,(2) 

AaA^j 

r)(3) 

Aa A^ 

r,(4) 

Aa Af, 

n(5) 

Aa A^) 

r,(6) 

yOyO 

Ad A^ 

XaX^ 

r,(8) 

AaA^ 



^1* ^2* ^2 ^1 I ^2* ^1* ^1 n2 

yairri/ ajkr^bikr^bjm ' r-^aim,r-^ajkr-^bikr^bjm 

— ( B^* (f'* if' + if'* B^* B'^ B^ \ 
yr-^aimr-^ajkr-^bikr-^bjm ' r^aimr^ajkr^bikr^bjm j 

ol* o2* o2 o2 I o2* ol* ol ol 
r^aimr-^ajkr-'bikr-'bjm ' r^aimr^ajkr'bikr-'bjm 

__ ( nl* n2* ol nl , n2* nl* n2 n2 \ 

y- aimi- ajkl- biki- bjm ' I- aimi'^ajki-^biki-^bjm J 

— ( B^* B^* B^ B'^ + B'^* B'^* B'^ B^ \ 
yr-^aimr-^ajkr-^bikr-^bjm ' r^aimr^ajkr^bikr^bjm j 

ol* ol* o2 ol I o2* o2* ol o2 
r^aimr^ajkr^bikr^bjm ' r^aimr^ajkr^bikr^bjm 

— ( B^* B^* B^ B^ A- B"^* B"^* B"^ B"^ 

\y aiml-^ ajkl-^bikrbjm ' l-^aimrajkrbikrbj 



nl* nU n2 n2 , n2* o2* ol ol 
r^aiml-^ajkr^bikl-^bim ' r^aimr^ajkr^bikr^bjm ' 
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APPENDIX E: F {iT Iji^ki^m) 
1. Matrix Elements 

The decay to neutrinos is mediated by charginos Xa ^ boson, and so the matrix 

element is 

2 

M = (El) 

a=l 

where 



\Pi Pk ) 

- ^ (7:..7a,™) u ipk) ^ PLt^ (p^) (E2) 

[Pi-Pk) -mi 

^ COS (pj — pj) — rriz ^ 



2. Total Width 

The width is therefore 



r (£r ^ IrvkVr^ = — ^ / \Mf 

JPS3 \ 6=1 



MzM%- 



r^^-;^- + 8km ^zz + 21^^- ) , (E4) 



where the partial widths are defined below. 

3- X~'X~ Contribution 

The chargino-mediated width is 



2 /-i-r? x^[\-x-rV'^ 



= E lliklajmlbiklljm dx X [l - x) (l - X - rfj f {x, Ta, n) 



TO,: 



a,b=l 
2 



5127^3 E lliklajralhiklljm^^-xl " (^5) 
0,6=1 

Here the integral ly-y- is the same as I'^o^a for the charge-preserving case, except the 
neutralino mass is replaced by the chargino mass and lepton masses ruk and are set to 
zero. 
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4. ZZ Contribution 



The Z-mediated width is 
„ m 



ZZ 



2567r3 



/ dx / dyMzMl 

Jx— Jv— 




2{x-l) i(') {x, rz) + 2 (l - X - rfj i^^^ (x, rz) 



(E6) 



Here again the integral is the same as the ZZ integral for the charge-preserving case except 
the lepton masses ruk and are set to zero. 



5. Zx Contribution 

The Z-chargino interference term is 



2567r3 



dx 



dyMzM 



5127r^ cos 9w 



Re 



I dx 2{l-x)g^'^'> {x,rz,n) + 2{x-l + r'fjg''^'> {x,rz,n) 

J X — 



rrii 



5127r3 



Re 



cos 9 



w 



Ihiklb 



ijm 



(E7) 



Here the integral Izx' same as J^^^o for the charge-preserving case, except the neu- 

tralino mass is replaced by the chargino mass and lepton masses ruk and rrim, are set to 
zero. 



APPENDIX F: P {ir ijqkQrr) 
1. Matrix Elements 



The decay modes with daughter quarks are much the same as those for daughter leptons 
except the neutralino intermediary contribution is removed and the different couplings are 
substituted. The matrix elements for the up-type quarks are 



Mz 

Mho 



(2) 



2 cos 9w {pi - Pj) - m 
-?,iaf^y^^^ cos9h 



u 



{Pk) {f, + (Pl - ^ sin^ 9w^ V {pm) Skm (Fl) 



V2 



{Pi-Pjf - m 



-U {pk) V {pm) Skm 



(F2) 



hO 
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M 



HO 



Q- (•'^1 (") • a 



Mao = 



Mgo 



-^i'^fivt^ cos/3 



{Pi-Pj) -mjjo 



U {Pk) V {Pm) Skm 



V2 



(Pi-Pj) -m\o 
1 



U {Pk) l^V {pm) Skm 



U (Pk) 7 V (Pm) 5km , 



a/2 {pi-pjf-m% 
where the factor of 3 is the color factor, and those for the down-type quarks are 



(F3) 
(F4) 
(F5) 



Mz = 

Mho = 

Mho = 

Mao = 

Mgo = 



1 



2 cos 6^1^ {pi-PjY 



u (Pk) [f- + pj) sin 9w - Pl)v (pm) Skm (F6) 



\/2 {Pi-PjY 



m 



-U ipk) V ipm) Skm 



hO 



-^if^jfVk^ cosOh 



1 



V2 



(Pi-Pj) -mjjo 
1 



U (Pk) V {Pm) Skm 



V2 {pi-PjY 



m 



U (Pk) 7 V (Pm) Skm 



AO 



-^i'^fiVk^ cos/3 



-U {pk) l^V {pm) S, 



V2 {Pi-Pjf-ml 
The total decay width is then 

r (i^r ijqkqk) = ^zz + ^hh + '^'^zh- 

2. Up- Type Quarks 

The decay width to up-type quarks is 



km 



(F7) 
(F8) 
(F9) 
(FIO) 

(Fll) 



r [i- ^ i-UkUk) 

9 



drrii 
5127r3 



9(t 



(2) 



COS 9 



w 



zz 



t=l 



"^ 2567r^m- ^^'^°^°^^>'°f^° "^horoIhoho + '^^hohoIIpho 

-\--^AOAoPaoao + ^GOGoIgOqo + 2A^ogo/^ogo) 



+2 



{^ZhoIzhO + ^ZRoIzRO + ^ZAoI^zAO + ^ZGO^zGo) 



COS ^^VK 

with integrals from Appendix C with quark masses substituted. The coefficients are 



, (F12) 



A 



(1) 



32 



- — sin ^w- ^ sin 9w + I 

/l(2) _ 32 4 8 2 /J 



(F13) 
(F14) 



22 



AhORO 

Aaoao 

Agogo 
Aaoqo 

Azho 

AzHO 
AzAO 
AzGO 



(2) (w) o 

(^ji Vk COS Oh 

(3) (n) . n 

tie CTjj (jjj sm t/^i^ cos 

(4) (u) Q 

ay vl' COS (5 



(5) («) ■ o 

Re 



(4) (5)* {u)2 . n a 

(^ji (^ji Vk sin/5cos/? 



= -Re 

= -Re 

= -Re 

= -Re 



— ^ 



(2) (w) n (° ■ 2a 

ayvl 'cos Oh ( - sm 9w 
ayyl sm 9h ( - sm % - 1 

" (4) (u) Q 

ay yl' cos (5 

(5) {u) . n 



(F15) 

(F16) 
(F17) 
(F18) 
(F19) 
(F20) 

(F21) 

(F22) 

(F23) 

(F24) 



3. Down- Type Quarks 

The decay width to down-type quarks is 

9m. 



r [ir ^ ijdkdk) 

9 



5127r3 



COS ^ 



w 



/ . ^zz^zz 

t=l 



+2 



(2) 



/TT 7— {BzhO^ZhO + BzHoIzHO + BzAoIzAO + BzG^^zGo) 

256v27r'^ cos 9w 



where again, the integrals are equivalent to those defined in Appendix C with quark 
substituted. The coefficients are 



B 



B 



(1) 
zz 

(2) 
'ZZ 

BhPh? 
Bhoro 

BfiOffO 

Baoao 
Bgogo 
Baog" 



gsm^^ 3 
4 



^ sin^ 6'w + 1 



9 



sin^ 9w sin^ 9w 



'^ji Vk sm^H 

(3) id) a 
(^ji Vk COS 9 H 



-Re 



^(2)^(3)* ((i)2 .^ a ■ 



(4) {d) ■ n 

ayyl'smp 

(5) (d) a 
ay yl' COS ^ 

-Re 



(4) (5)* {d)2 . r) o 

Vk sm/3cos/3 



(F25) 
masses 

(F26) 
(F27) 
(F28) 
(F29) 
(F30) 
(F31) 
(F32) 
(F33) 
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Bzh° 

BzH° 
BzAO 
BzGO 



-Re 

Re 

-Re 
Re 



aj-^l/J,'^^ sin Oh sin^ dw - 
of^yf' cos Oh(^ sin^ - 1 

" (4) (d) • n 

(5) {d) n 

ay yl' COS 13 



(F34) 

(F35) 

(F36) 
(F37) 
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